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NOTATION 



m 

A, B 



C i 



e 



e rr* e zz* e ee 

E 



G 

L 

r 

T 

u 

U 

V 



w 



' rz 



°zz , 0 rr’°ee 

T rz 



maximum rigid body acceleration of the projectile 
material constants, eq. (8) 
constants due to body forces 
dilatation 

normal strain components 
Young's modulus of elasticity 
body force in the z direction 

coefficients of eqs. (21) from Rayleigh-Ritz method 

shear modulus 

length of propellant 

radius, inner radius, outer radius 

total potential energy 

displacement along r axis 

elastic strain energy 

potential energy of external forces 

displacement along z axis 

Rayleigh Ritz undetermined coefficients 

shear strain component 

weight per unit volume 

normal stress components 

shear stress component 



INTRODUCTION 



The launch phase of a gun-launched rocket is associated with high 
acceleration loading of up to 10,000 g. As a result of the high accel- 
eration the motor case portion of the rocket, which contains the 
propellant, is subjected to axial and lateral loads. The lateral loads 
are associated with the deformation of the contained propellant. This 
"slumping propellant" behavior is the concern of the present investi- 
gation. 
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ANALYSIS 



The propellant-motor case system may be idealized as a two-dimensional 
system because of axial symmetry. With cylindrical coordinates, r and z, 
a typical section is shown below. 



fore 



z, w 



aft 




u = displacement in the r direction 
w = displacement in the z direction 
r^= inside radius of propellant 
r Q = outside radius of propellant 
L = length of propellant 
The particular problem of the unbonded 
propellant is considered here. The 
stiff motor case shell is assumed to 
completely restrain axial displacement 
at the aft end (r,0) as well as radial 
displacement along the outer radius 



(r 0 ,z) with the resultant boundary conditions. 



w(r, 0) = 0 
u(r 0 ,z) = 0 



T rz (r,0) = 0 



T rz^ r n» z ) = ^ 



The remaining sides (r, L) and (r^,z) are stress free, that is 
o z (r,L) = 0 T rz (r,L) = 0 
o^.z) = 0 T rz (r i ,z) = 0 



(la) 

(lb) 

(lc) 

(l d) 



The acceleration of the propellant subjects every particle to the body 
force 

f z - -pa m (2) 

where p is the weight/unit volume, and a is the maximum acceleration 

m 

in g's. 



2 



A solution to the problem is obtained by the Rayleigh-Ritz method. 
According to this method, an approximate solution satisfying the displace- 
ment boundary conditions is formed. The unknown coefficients of the solution 
are determined by minimizing the total potential energy of the system. 

Here the displacement fields 

u(r,z) = (L-z ) (a + er + yr 2 + 6r 3 ) = (L-z)f(r) (3) 

w(r,z) = z(L - |-)(n + \r + £r 2 + cr 3 ) = z(L - |-)g(r) 



are assumed. The linear z term in the u displacement and the quadratic z 
term in the w displacement were chosen to satisfy the following conditions: 
i) Results from the Rohm and Haas finite element solution, 

ref. 1, Appendix B, for the accelerated, unbonded propellant 
show a r and o z to be linear functions of z. Note that Equations (3) 
give this condition, i.e. 



o 



r 



a 



z 



E | 


3u 




IF 


E 


3w 


1-v* 


3Z 



L 




E 

Vv 7 




(L-z) 



E 





ii) The propellant is assumed to be nearly incompressible, and thus 

linearly varying compression along the z coordinate causes the pro- 
pellant to move towards the centerline since it is restrained by the 
case from moving outward, i.e. u $ 0 for all z,and u = 0 for z = L. 

This assumed displacement field is not valid for the bonded propellant. 

The equation for w satisfies the displacement boundary condition w(r,0) = 0. 

The displacement boundary condition, u(r Q ,z) = 0 is satisfied by setting 

a = -Br 0 - yr 0 - 6r Q (4) 

The total potential energy of the system, T, resulting from the displace- 
ment fields given by equations (3) must be formed. Denoting the elastic strain 
energy as U and the potential energy of external forces as V, we have 

T = U + V (5) 
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Determination of the elastic strain energy U proceeds as follows. 



According to the strain-displacement relations 

e r = fr = ( L " + 2yr + 36r ^ 

e e “ r * ( L ' z ^r + e + yr + 6r ^ 
e_ « |S- * (L - z)(n + xr + + ?r 3 ) 

Z d Z 

3 U . 3 W t , , „ 2 . . 3 \ 

T rz 3z 3r ' 1 

I 2 

+z(^r - z)(x + 2^r + 3cr ) 

and the stress-strain relations, 

°f * fle r + B(e z + e 0> 

°r = Ae e * B(e r + e z ) 

o 2 = Ae z + B(e r + e 0 ) 

t = Gy 

rz T rz 

where 

A = frf jjrf— 

R - vE 

B - (1 + v)(V - 2v) 

and 

G = 2Tr+"v'r 



( 6 ) 



(7) 



(8) 



(9) 
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Substituting equations (6) into (7) yields stress-displacement 
equations. The elastic strain energy of the system is given by 



u 



°u dV = i 






(o r e r * °e e e + °z e z * T rzVz> dV (10) 



where the integration is over the volume of the propellant, 
independence with the e coordinate gives 



The 



U = 2rr 



. 1 



r i *- 



r o 0 



(o r e r + °e e e + °z e z + T rz Y rz } r dr dz 



For convenience, let 



(ID 






o e r dr dz 
r r 



U b = 77 



o 0 e 0 r dr dz 



( 12 ) 



u c = » 



o z e 2 r dr dz 



u d * * 



VzVz r dr dz 



then 



U = U + U. + + U . 

abed 



(13) 



Substituting the stress-displacement, and strain displacement relations 
into equation (12), and integrating, we obtain 
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II - 

u a - T" 






A {?A x+ i TV+ i r/l 









4B>{o<y9/V 4^* (/6i| -*-y6 i + a.To l) + ( 3 yS 2T 

+ a1\43*5«c +-yS\) -»- ik 4 (,\p& 4^+34^ 

+ A?X 4^0 4 (5^ + 35X4/9^ 

O ' — « ' 

+ **0 +-f- 6 (3<5 i -V-3i^+5LYc)- 1 .^ < S?A 7 | 



(14) 



Uu = 



3 

TfL 

“ 



A «>£ Z X*v/l 4 3>«</SA, 4 ~ ( |3> X i 3ot^ 

+. A 3 ( a°t 5 4 Zj3 if) + ^ ( t X 4 Ay8 <s) 

+ |.y«5/t 5 4 ±<5 x a‘J 

•f B ^ ( 0<Y l + <X ) 3 ) ^ (a.ot8'+jQr^4y3 :> '+o(X) 

+ iL B ( 3 ot^ 43£tf 4*^ 4/3X 4 - 0 * 5 ') 

+ 4. 4 (4^3<^4 ajrSSti + TX 4/8f 4*0 
+ 4? (stfS 4 5 X 4 y ^ 4y0^') + f t (3i :L -i- < S^4Tr^) 

5 7l ft 0 

*“*L 



u = 



- ^L 1- 



A\-i ^4 -J^X n*4 £- 4 (X* 4=^)4 ^?(a.AS + a ‘ 1 l^) 

^(^aA^l + i KA 7 + <st|! ^ 

4 B ^<**1^ + ^ (^y 9 ^ 4 04 x) 4 ^(3 yy^4A./3>i + * 5 ) 

+ I 4 ( 4 ^ + 3T)i4a / 3^4o< ? ) + 4 s ( 4 <SX 4 3fff 

+ ^ 5 ) 4 ^ (4if + au)+ 1 A 7 y| * 

-*AL 



(15) 



(16) 
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U d = * G + f ^ C/6 a, + ao<r)+ 

+ 2.Ty8) + ^£<5/1 ^ + y TJNjA 7 } 

(*<*§ + £*) + ^ (3«»<C +^5 
-t-^A) +- -|? (3y8c +AV§ 4- 5A^ + l fc ('3t'C 
+ Af<0 +|- i^A 7 } 

+ roU xV% T x§A3 + A " ( ^ + i A ^ + T^ A ' o (17) 

+ i* v }| 



The potential energy of external forces V, associated with the inertia 
body force f given by equation (2) is. 



V 



r o L 



f • w dV = - 2. 



p a„ ’ w ’ r dr dz 
m 



r i 0 



(18) 



Substituting the expression for w from equation (3) into equation 
(18) yields 



irL - 2 2 3 l l *2 5 -, 0 

v = ‘ T - p a m ' [nr + I Xr *7 5r + I ;r 3 

r 1 

The total potential energy of the system is obtained as 



(19) 



T = U + U K + U + U . + V 
abed 



(20) 
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The theorem of minimum potential energy states that equilibrium is 
associated with a minimum of T. Here T is a function of the 
independent parameters B,y,<5,n>S and c. The minimum of T is obtained by 
taking the partial derivative of T with respect to each of these seven 
coefficients. This yields the system of 7 linear algebraic equations 
in 7 unknowns. 
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The explicit expressions for the F. . and C. (i.. = 1, 7) coefficients 

1 J 1 1 J 

are 

F(',0 - A {-ZAi-l/il 4-^/1^; + Mi 1 

+ B {- } + 3 G { i. n?Al 4- ±*0*1 -1a!} 

= Aj-71 0 3 -a/l; +/) 0 Jll J u ^0 J 

+ 1 - 1 / 1 , A* } + 3G Xf z flS-V Jii 

, 3 l 30 

/1 0 4- A 0 *x - J } 

F(V) - F(‘A) 
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p(i,3)-fU 0= A lr a ^"^t + +a/i 0 4 JU^| 

+ B ^a 0 a a 3 - a +-*/i 0 3 4;} 

+ A^“ ^ o A* + J A < + J^o-^a " 3 *jL } 

f(i, 4) = F(4 y l) = B {-Ail; + ia^.I 



F(l,s)' F(5,0 - 



B{±.(a*-A;)-A« + \ 

-a* W^+i^-i**} 



F(.\,<o)~ F(fc / l)= B|^-aJ + jA 0 >i;-^A } 

-afc{~£/C + |A 0 ji- } 



F(l J 7) = F(7,l)^ 



B{i*MM!-fA!] 

- i6 l-|b /, » + l ;, ' J) *-T^l 



F(a,i) = A {-A, 4 - § At + *4 + a/l. 4 ^ | 

+ B^-i^+4< 

+ 3 G^/ { 3* Aq ~ A 'k + ^0 A a. “ "3 } 

F(a,a)^ F(3,0 = A{g/l* - ~ n 5 i + - 1 - n * A ^ Xvc ^T°) 

+ B {-**1 + **l A* + a.A>i' A M 

* %S \% n\- *» A? + W / ) J + 1 «J A? - I < } 

F(vA = F(4,l) = B{-h’w<A;| 



F(a,s) = F (S,a.) = 






f(i>1 = f(m-)= b { il ^ + j n i } 

F(a,7)“ F(7,l) = 



F ( 3,3) 




f(3>) = F(4,3)= B {- ini n.^ } 



F(3,5) = F(s, 3) = B{|-(/l, 5 -yi!)-4 t S +/l, 3 4j} 

- zG \- f o a! + i <1 3 j 

F( 3 ,«) = FU,3) = b{| y •*; ) 

F(3,7) =F(7,3) = B{-£.A. 7 + i*!*;- fJll } 

4*1} 

F(4,4)= A {n^-A? i 
F ( 4 5 ) = F(S,4)= A | -j ( <1^ ~-4i )} 

F(4>) =F(t,4)= * {{ U?~4;)} 
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F(4,7) = F(7,4)= A {a (-I, 4 -A; >} 

f(s,s) = + GL^f (<-a* )} 

F(s,t) = F(fc,s) = +GL l i-|(-1„ 3 -.1 x 1 )} 

F (s/7) = FUs) = h{%Ui- aJ)] + GL l {§• ( 'K'-A? )} 

F(<i>)- A{i(A‘-A.hJ + 6L !l {±.(^-Al)} 

F(t,7) = F(Z6)» + 

F ( 7^ 7 ) = A{±('>®-A*)} + 6L 1 {|- Ui-At)} 

C,= C»- C 3 =0 

4 ~ ^ A, - ) C 5 = - a * )/a . 

C 6 = C 7 =-ip<wU*-A')/s. 

A computer program was used to obtain the solution of the problem. A 
large number of problems were solved for various values of the physical 
parameters, r^ , r , L and a m , and the propellant properties v and 
E. A listing of the program is given in Appendix A. 

DISCUSSION OF RESULTS 

It has been noted by other investigators that numerical difficulties 
arise in the analysis of propellants when Poisson's ratio approaches 0.5 
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(ref. 2). To circumvent this difficulty Hermann (ref. 3) gives a 
special variational theorem for nearly incompressible materials. In 
the present investigation the numerical problem was overcome by using 
a CDC 6600 computer with double precision. This results in approxi- 
mately 32 significant digits, sufficient for Poisson's ratio to 
0.4999999. 

A study was conducted on the effect of Poisson's ratio on the 
propellant stress for a propellant with 1.5 inch outer radius, 0.5 
inch inner radius, 9.5 inch length, 500,000 psi Young's modulus and a m 
equal to 8000 g. Poisson's ratio was varied between 0.45 and 0.4999999. 
The computer program results for a r at the propellant-Case interface at 
the base are given in Table 1. 

Table 1 
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Figures 1 through 3 are curves showing these results; the first 
figure being a direct plot of stress o f versus Poisson's ratio v. 
Figures 2 and 3 show more clearly that the stress approaches a finite 
value of -4650 psi as Poisson's ratio approaches 0.5. In Figure 4 the 
dilatation e = e Q + e r + e z> at r = r Q and z = 0 is plotted against 
Poisson's ratio. The resulting curve verifies that the incompres- 
sibility condition e approaches zero as v approaches 1/2 is satisfied. 

In addition to the study on the effect of Poisson's ratio on 
propellant behavior, analyses were made to determine the effect of 
other parameters to behavior. The results obtained show that the 
radial stress varies linearly with respect to acceleration, radii, 
length and Young's modulus of elasticity. 

The results of this analysis were compared with two other sources; 
a Rohm and Haas finite element analysis (ref. 3), and an exact 
solution for an incompressible material with zero interior radius. 

The present formulation resulted in a 20 percent higher radial stress 
than the finite element program, and 2 percent higher than the exact 
solution, at the propellant-case interface at the base. 

CONCLUSIONS 

A computer program was developed for the determination of the 
stress distribution due to a slumping propellant contained within a 
launched rigid motor case. The program was initially incorporated 
into the SATANS finite difference program for the calculation of the 
elastic static buckling load of the motor case subjected to high 
acceleration loading. The stiffening of the motor case, which is 
predicted by theory, was obtained in the SATANS analysis; a favorable 
structural response. 
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An unfavorable effect on structural integrity is associated with 
the additional contribution to yield from the propellant radial stress. 
For a motor case of 2.5 inch outer radius, zero inner radius, 1/8 inch 
thickness, 20 inch length, and a maximum acceleration of 8,000 g, 
the radial stress of approximately 10,000 psi gives a hoop stress of 
100,000 psi. The contribution to yield of this stress component is 
of significant magnitude. That is, the propellant has a beneficial 
buckling stiffening effect but an adverse affect on yield. 

The present analysis is based on the linear elastic theory of 
solids and, therefore, the results are not valid beyond the initial 
yield point of the propellant. Because it is difficult to obtain the 
material properties of propellants (Young's modulus, Poisson's ratio, 
and initial yield), definite assertions regarding the structural 
integrity of the propellant are not given here; however, the stress 
results obtained in some of the analyses show that additional study 
regarding this concern may be warranted. 
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